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T'HS UoE OF j;iB;.OUTl!;SS V.ITH TiiE BENDIX G-l5 



Tnis paper is •j.nteiidsd ;'S a gtiide for follcying the Specification Sheets 
of tlie subroutines published for use v^ith the Bendix G-lS. 



Iti 



its execution. Caj/e s^/ould be 'taken tl.c 
must be stored (te-apora.rlly, at least) S 
to be executed V 



:..:.::.■. •■^..he subroutine rriust te stored for 

h v.^;^:..-; '-i subroutine is tc he used, it 
a the co;ruuarid line fro^u vidch it is 



2 . Input ? 

a. Most su:..routi..nes specify thai tlv- :input to the subroutine^ sucu as an 
argument (for instanc-, an an^le of vrhjc/i a trigonometric function is desired) 
ce placed in a pre-deteri^vined stora^-e locatj.on, such as a Kord or words of the 
last-access lines or \.:.o-word regi^tars, lefore transfering to the si±^routine. 

b. A '^lietuTu Cornnand"' ma:/ be required iri the AR on entering the subroutine. 
T^-rls ccnn-aand v;ill ce >:>;:-cuted 'ij)i the subroutine when its fimction is cOT^vpleted, 
to return control to '.he main routine- \Zvb ¥Bxk. and Ret'urn featiare of changes 
of -onuKand line in '.he n*l5D are not ':.sv;d in subroutines published for tliat 
riodel^ All Return Cv.nrands should be written as l^arked I^^ansfers,) 

h. En tr:/:; 



This refers to the word locatio:. in the suc>routijrae cormnand line to which 
control should be ".x^ansfered. Give a con;^:'^and to transfer control to t!ie 
co--mrid line of th.-^ s iil^routi^ne , with tliis location as its ^'N"*. (For G-ljD, 
a har;-.d Transf-r. } 
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U^ Exit : 

The Specification Sheets will give tlie word tjjne at which the Return 
Command will be executed* The "N" of the Ret\irn Command should specify the 
location and (for G-l^A, the source) (for the G-l^D, the characteristic) the 
command line of the main routijie to wM.ch control is to be returned. 

5* Scaling g 

Nuiabers In the G-l5 are generally scaled so that they appear as fractions 
(less than l.O) within the macliine. Some subroutJJies may specify whether 
their inputs (and/or outputs) are to be scaled with either a binary or 
decinal scaling factor* 

6. Short Line Use ? 

ifeny of tlie Subroutines will make use of the fast access lines* Generally 
such use will be noted in the Specification Sheets. Informatn^on stored there 
by the programmer prior to subroutine entry will have been destroyed by the 
execution of the subroutine. Tlfie information if retention is desired, should 
be placed in storage locations not used by the subroutine. 
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THE USE OF SUBROUTINES WITH THE BENDIX G-l5 
TITLE 




3AHPLE PROGRAM — USE OF SUBROUTINE 



PaOBLE>h To find tiie natural log of x. 

GIVEN: Log subroutine in Line 8 

X less than 100 j x • 10"2 in 20.02 



Word Position 
of Command 



00 



06 



EXPUNATION 

Loading the odd side of ID from an even location, 20,02, takes 
two word tJnes, 02 and 03, and is done double precision. This 
clears the even side of ID and all of PN. x is at a decinial 



-2 



scaling factor of 10 



2 «7 

The conversion factor, 10 • 2 S is transferred into the !^ 

during word position 0?« 



09 



2] 



26 



Since all binary digits of the multiplier after the fifth are 
zero, only the first five bits need processing, requiring 
nultj.pli cation for ten word tiiaes. Multiplication by the 
conversion factor has eliminated tiie decdmal scaling factor 
and introduced a binary scaling factor, 2"^, leaving x in a 
form suitable for ttie log subroutJ-ne, 

x»-10'2^ vtnch is needed for later calculations, is moved from 
20.02, where it would be destroyed by the subroutjjne, to 21.02. 

i^'2^^- is moved fi'oia the odd side of PN, via AR, to an even word, 
20.00, Since an odd v7ord tdiiBS, 27, and an even word time, 28, 
which are not parts of the same double precision word are jjri- 
volx^ed, the command j.s written as a single precision block 
command. 
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Word Position 
of Command 

29 



36 



58 



u6 



03 



EXPUMTION 

A block command for 108 word times transfers the Log Subroutine 
from Line 8 into line 1 for execution. 

The subroutine requj.res that the exponent of the binary scaling 

factor of X be placed in 20.01 as an additional input. This 

-28 
exponent is itself expressed at a binary scaling factor of 2 

The "Return Command" is placed in the AR prior to entering the 
sxabroutine. It will be executed from LOl.Ol as an exit from 
the subroutine. Since its "N" is 03 > control will be returned 
to LOO • 03 in the main program • 

L01»00 is the entry point for the Log Subroutine, single pre- 
cision* The "N" of this command is 00. 



Log X at a binary scaling factor of 2^^ is placed in the AR, 
from which it may be typed or used for fvtrther calculation. 
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IT 



N«. 1101 



DESCRIPTION 



This subroutine consists of a Binary to Decinal subroutine, a Decimal 
to Binary subroutine an aji output subroutine. The first tvo are 
standard conversion routines with normal entries and exLts, as described 
on their specification sheets. 

The output routine scales the nunb^r then converts It and types. If 
the scaling is perfprmed by multiplication, the scale ftictor (s.f.) 
must be stored vith a negative sign, e.gt - (s.f.) in cell 21.00. 
If the scaling is by division, (s.f.) is stored in 21.00. This 
in no vay affects the sign of the number being processed. If a 
new format is to be used, it is loaded into IDq i cuid the routine 
is entered at word time uO or u6. Othearwise entry is at word time 
05 or 06. Typeout may be executed from either line 19 or from AR 
as shown on specification sheet. Line 19 typeout allows computation 
not using line 19 to proceed. 
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SPECIFICATIONS 

l^e SubroutJne (G-l^D) 

Mode of Operation Plxed Point, Single Preci sion 

Execution Line 02 

Entry !^ypeout 



!^ypeout 
from 


Unchanged 
format 


New 
format* 


L19 


05 


u6 


AR 


06 


uO 



Exit At^word time 37 

Data Input t \'^\ — >- 21.01 

rtet. Coram. ->-AR 



Scale Factor 



-^ 21.00 



♦Format 



■^ID, 



'0,1 



Data Output Typeout under control of format from 

line 19 or from AR, 

Executdon Tjme. . Line 19: 3 ftev» to initiation of typing 

Computation may proceed while typing* 
Afti !iypjng completed before exjt from 
subroutine. 

Error Stops Bell rings jf improper scale factor used. 

Short Lines Used 21.00 and 21.01 only 

Check Sum Cell Used u? 



Tl^/70^y 



Compi/ft 



TWPUT-QUTPU T ^^'l^^^^ '^ DKHTMAT,^ 
Prtparad by : Ed Williams 
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Pan; 1-30.^7 
N«. 1101 



SPECIFICATIONS 

lype Subroutine (G-l^D) 

Mode of Operation FJxed Point, Single PrecJsion 

Execution Command Lj.ne 02 

Entry Word time 61 

Exj t Vford time 63 

ScalJng Convert x as a fraction 

Input X (binary) in ID-^^ 

Return command in AR 

Output 1 1^1 i (decimal) in AR (7 djgits and sign) 

Execution Time 1 revolution 

Error Alarms None 

Short Lines Used None 

Check Sum Cell Used u? 

REMARKS 

This routine, the Decimal to Binary and the Output routine are together 
as a una tin the same 1 J ne . 
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SPECIFICATIONS 

Type Subroutdne (G-15D) 

Mode of Operation Fixed Point, Single Precision 

Execution Conanarid line 02 

Entry Word time U6 

Exit i , Word tdjne U? 

Scaling FractJ.on or Jnteger 

Input X (decimal) in ID^^ (7 digj ts and sign) 

Return command in AH 
Output 1 1^1 1 (binary fraction) in MQq 

+ |x| t (binary integer) in AR 

Executjon TJme 5 revolutdons 

Error Alarms None 

Short Lines Used None 

Check Sum Cell Used u? 

REMARKS 

Thds routine, the Binary to Decimal and the Output routine are together 
as a unit in the same line. 
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Format 

Clear Line 19 
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Format 



Set Switch 
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Set Switch 
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Set Return Link 
Scale Factor — >-AR 
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X by |B.f.| 
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X' to Deciinai on page T 
r xa 
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Exit 
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* All Commands are on page 6, with this exception ^nly. 
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TITLE 



?N CORPORATION 



Prenr«d by; H. LeTds 



LOS ANGELES 45. CALIFORNIA AP>fO»»l by '._ 



■f 



..^^^nt^rv 



€>'^^*'yL^<^X^ 



Pail 1 of 15 
Date: 3-13-57 
.No. 1202 



SPECIFICATIONS 



Types Subroutine for G-.15D 

Mode of Operations • • • • • Fixed point, single precision 

Executions ••<» From Command Line 1 

Entry? « * At word time 093 

Exit ooooe^o ^Q't* ^0T^™« ^rom Line 1 at word time 098 

Scaling? «•.••.•»•••. ♦••B must be in circles; 9 circles m 

6 degrees = 6 radians , where -1 <^ 6 <^1 
3S0^ Ttt 

Data Inputs . . « o © >-20,02 

Ret, C omm* ^ AR 

Data Output? , l/2 sin 9 r 20.00 

e . 20.02 

Ret. Comm, . Loc. 098 of Line 1 
20.01 is undisturbed 

Execution Timet . o 203 msec (7 rev.) 

Error Stopss ao..oo«««c««o None 

REMARKS 

Command in location 050 can be changed to clear (Source 23, Destination 31) 
to prevent obtaining an answer -0, 

All of the short lines are used in this subroutine • 



SINE 



(G-15D) 



TITLE 



Date : 

5-13-57/ 



Pate 2 ot 05 
No. 1202 



METHOD 



The argvunent, 9 , must be in circles , 

e circles a 9 degrees n 9 radians 

TIT" 



"sSSS" 



Consider © circles - + 



K -4- 

— V 
4 -^ 



, where K is the number of right 



angles mod, 4, and y is the acute angle where ^ y <^ ^, 



2 



Jl 9 



, we use the same subroutine for 



Since cos 9 - sin 
both functions. 

The following table shows the instructions used to develop the proper 
sign of X and the proper X for the polynomial evRluation, 

5 

b 



ir sin 



V 



TT X - X 2 

2 iso 



-X 



2\ i 



where: 



-1<X< 



1 and where 



■ b, are the positive constants. 
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TITLE 

Prepared by : H» Lewis 



LOS ANGELES 45. CALIFORNIA ApprOVed by '■ / M^^>t-^>^^ 



^ 




SPECIFICATIONS 



Typet •••• Subroutine for G-15D 

Mode of Operation; Fixed point, single precision 

Executions •••««••• From Command Line 1 

Entrys .At word time 097 

Exits <,**•* Ret, Comm. from Line 1 at word time 098 

Scaling^ 6 must be in circles; 6 circles m 

e degrees » 6 radians , where -1<C 6\1 
5g^5tr" 2TT~ 

Data Inputs « - © ^20.02 

Ret . Comm. ^ AH 

Data Outputs l/2 cos © r 20.01 

e m 20.02 
20.00 is undisturbed 
Ret. Comm. s Loo, 098 of Line 1 

Execution Times 203 msec (7reY.) 

Error Stopss o None 

REMARKS 

Command in location 050 can be changed to clear (Source 23, Destination 31) 
to prevent obtaining an answer -0. 

All of the short lines are used in this subroutine • 



COSINE 
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TITLE 



Date : 
3-13-57 



Page lx_ of _15 
No. 1202 



METHOD 



The argximent, © , must be in ciroles. 

e circles - © degrees i 
5555— 



e radians 
^__. 



Consider 6 oiroles r 4- 



£+ y 

4 



, where K is the number of right 



angles mod, 4, and y is the aciTFe angle where <^ y<^4* 

Since cos 6 r sin( "TT + q \ ^ ^q use the same subroutine for 
both functions. 



The following table shows the instructions used to develop the proper 
sign of X and the proper X for the polynomial evaluation* 



i sin 



TT X « X 



t s 



-X 



where: -1<C X<f 1 and where 



ai 



hi are the positive constants. 
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NOTES: 
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STEP 7 IF NOT SET, MEANS ON y-AXIS FOR 


SIN e, X-AXIS FOR COS e. 
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SIN 0, y- 


AXIS FDR COS e. 
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FORM 
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SINE, COSINE 



(G-15D) 



TITLE 



Date: 
3-13^7 



Paje 2__ of 1^ 
No. 1202 



FLOW CliilKT 



ENTER SIN SR 
AT Ol..-.- 



RET. COMM. 



ENTER COS SR 
AT '.-l.yl 



RET. com.- 



01.98 



01.03 



SET UP STORE COMM. 
(^ sin e , PN-, ^^^. 20.00) 



01.09 



SET UP STORE COMM. 
(^ cos 9 . PNn ^20.01) 



01. u3 



SET AR - 



'' 01. ul 



SET AR ■ i 



01.24 



DETERICNE SIGN x, THEN 2©" 
Ref. Specifications 



RESET OVERFLCW 



C 



r 01.45 



29" ^TEST 



C 



> 



-0- 



01.47 



2 



-PNn 



jf.O 

V 01.48 



01.50 



26"- 



-TEST 



> 



,0— ^PN^ 



I 



01.51 



DETEmnNE ±x 
Ref. Specifications 



01.53 



LOOP TO EVALUATE 
POLYNOMIAL 



- Z \ (-') 



2,1 ■•■ 



isO 



01.92 



EXECUTE STORE COM. 



EXECUTE RET. COMM. 




SQUARE ROOT 



(G-15D) 



SSDBD 



masmm 



Rffl 



TITLE 

Prepared by : h, Lewis 



LOS ANGELES 45. CALIFORNIA ApPTOVed bjf 



-^4^^ 



./ 



Date: 3-13-57 



No. 



1202 



SPECIFICATIONS 

Type? Subroutine for a«15D 

Mode of Operation: ••o*»,«# Fixed point/ single and double precision 

Execution? **,o««oo».«* From Command Line 1 

Entry? a«o«««. • At word time 094 

Exits •«o»»».0«*»«.» Ret, Comm, from line 1 at word time 098 

Scalings • o . , N< 1 ; W<C !♦ 

Data Input: o . . • * * N "*" > PNq -j_ 

Ret, Comm. >- AR 

Data Output: « ^"double precision « PN - ; 

1F"single precision ■ 20,63 
II s 21,00; 21.01 

Ret. Comm. r Loc. 098 of Line 1 

Execution Time: ,,,*,,,«• 261 msec (9 rev.) 

Error Stops: , » • ♦ • If -N (double precision) is entered, 

routine will ring bell and halt. 

REMARKS 



Programmer must test for -N (single precision) since routine does not detect this^ 
All of the short lines are used in this subroutttie. 



SQUARE ROOT 



(G-15 D) 



TITLE 



Date: 

13-13-57 



Page 9 of 15 [ 
No. . 1202 



IJETHOD 

The routine finds the double precision square root of the number, N, 
But it is written so that the arcsine and arccosine routines can use it. 
to find a single precision square root. So first, the single precision 
square root is found as follows: 

(l) r_^ a ■!■+ ■!• N is used as the first approximation, 

where N<r 1. r^ vdll be greater than the square root. 
Then ■'.'.'e use the formula, 
^^^ ^1 = N -*- Tq and iterate until ( r3_ - Tq -+• 2"^^ ) ^ 

Now that we have the single precision square root, two double 
precision divisions using formula (2) are sufficient to form 
the double precision sqxiare root. 



FORM 
BCD 

108-7.0 



SQUARE ROOT 



(G-15D) 



TITLE 



Da t e : 



Page lQ_of J5 
No. 1202 



FLOW CHART 



NOTE 

Sin"-"- & Cos-"''" 
Enter at 001 



01.58 



-*- 20.03 



NOTE 



> 



-1 -1\ 

Sin & Cos 

Exit at 

037 & 033 



ENTER SQ. RT. SR AT 094 
SET COMM. »-2l,03 



(PN, 



0,1: 



-TEST) 



01.96 



RET. COMM.- 



Lgs 



'■ 01.07 



(^TEST SIGN "i-J"^- 



i 



01.49 



{^"N" ^^-TSSt) 



^0 
t 01.59 



r. = *+ * N 



SET AR NEG. 



01.85 



r^ »-> 20.03 

SHIFT TO FORI,! r, /g 

\ 01.88 



(^TE3 



^ 



TEST SIGN AR / 

t "^ 01.92 



EXECUTE COim. IN 21.03 



'N" 

■iwdftii 



i. 



— ( 21.05 s PNq^i 



-TEST 



) 



01.28 



HALT AND 
RING BELL 



"N" :^ 01.99 



FORM 
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DOUBLE PRECISION 



^1- 2ro+-2- 



r - JL- H- ^0 



2 = 2ri 



*■ EXECUTE RET. COM. 



01.91 



SINGLE PRECISION 



2r<, 2 



SET 



ARq 



- r.- 



SET r Q s r j^ 
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TITLE 

Prepared by : H. Lewis 



LOS ANGELES 45. CALIFORNIA ApprOVCd fajf '. 



^ 

J^ 



k 



^ 



Ai-^ c^ 




Pan ja oi 15 

Diu: 3-13-57 
No. 1202 



SPECIFICATIONS 



q cf degrees 
= 5600 



<i radians 
2 ft 



1. 

4 



Type? ««*o«.»*««#«** Subroutine for G-15D 

Mode of Operationr •••»♦•# Fixed point, single precision 

Executions <,•»•* From Command Line 1 

Entry: . . • * At word time 078 

Exit: • • Ret. Comm. from AR at word time 054 

Scaling: . * . . ^ r sin" 

^ \ <sin-^ S< 
Data Input: \ S v-20.00 

R et » Comm, ^ AR 

Data Output: s S - 20*00 

\ cos ^ 2 20*01 

<^ » sin""^ S m 20,02 

Ret, Comm, s 21,02 

Execution Times • , * o -. • * * • 754 msec (26 rev*) 



Error Stops: ,.««6*.»«« If \ ^^ ■ ^^ 

I 2 



- (fi)2 



is negative, routine will 



ring bell and halt* If \\ ^'^\\ 
routine will ring bell and halt. 



\ -.(Ir )>i. 



RE^IARKS 



Tan" (y/x) subroutine must be stored in Line 04, After use, tan" SR is 

left in Line 1, 

All of the short lines are used in this subroutine. 



ARCSINE 



(G-ISD) 



TITLE 



Date: 
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No. 1202 



METHOD 



Since sin"l S - tan'^ 



J first 



\ 



8^2 



X- (-) is evaluated 



'Nl'i-s"^- 

by the square root routine, then assigning ^ S r y and 



i- (1)^ - X, 



■W9 enter the tan"^ ( - ) routine to find fi . 



FORM 
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TITLE 



Prepartd by : H« Lewis 



LOS ANGELES 45. CALIFORNIA ApPTOVgd 
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•^^^>7^<r-^/ 



/. 




Piia Jjol li 
Date: 3-13-b7 



No. 



1202 



SPECIFICATIONS 



cf degrees 



C^ radians 



Type: •♦• Subroutine for Q-ISD 

Mode of Operation: Fixed point, single precision 

Execution: . , • • • • From Command Line 1 

Entry; ♦ At word time 073 

Exit: Ret. Comm, from AE at word time 054 

Scaling: « ^ t cos" G s 

<^cos""^ c;< i- 

Deta Input: . ^C 3-20.01 

Ret, Comm. 3^.AR 

Dat)^ Output: i sin ^2^ ~ 20,00 

I C s 20,01 

m cos"*^ C r 20.02 

Ret, Comm. ^ 21.02 

Execution Time: 754 msec (26 rev,) 

,^ 2 

Error Stops: •••« If 



i - ( ° )' 



is negative, routine will 



ring bell and halt. If (|- C ;, 



±i- (If) >1, 

routine will ring bell and halt. 



REMARKS 

Tan""^ (yA) subroutine must be stored in Line 04. After use, tan SR is 
left in Line 1« 



All of the short lines are used in this subroutine. 



ARCCOSINE 
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TITLE 



Date : 

3-13-57 



Page jl^ol a^ 
No. 1202 



METHOD 



Since cos" C - tan 



-1 .^ll^ . first^Ji- (o/2)^ is evaluated 



by the square root routine, then assigning^^ (o/2) - y and -I* C s X 
we enter the tan"" ( ^ ) routine to find ^» 
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TITLE 
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SET a m 
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ENTER SQ. RT SR AT 001 
FOR Vn, single PREC. 



EXECUTE COMM, in 21.03 
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01.57 



SET 



y-A/i-(#) 



01.60 



ENTER TAN""^ ( I ) SR 



AT 060 



-1, 



FOR J2( - SIN'^S or COS" C 



EXECUTE RET COMM. 
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TITLE 
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LiFORNiA Approved by:. y^A^L^^,c^./(^nL 



^ 



Paje _l_ol J_ 
Date: 3-13-57 



No. 



120b 



SPECIFICATIQNS 

Typar Subroutine for G-l^B 

Mode of Operationr . , . . Fixed point, single precision 

Execution?* Command Line 1 

Entry t •#•#•••••• At word time 059 

Exitr * 6 . . ^Bet* Coram, from Line 1 at word time 067 

Scaling? The argument and all tabular values of x 

are at one scale* All tabular values of y 
are at one scale, not necessarily the same 
scale as the x values. 



Data Inputr • , x 
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-vAR 



Pickup Comm* (deferred, T « "t;"^ j^ - ^%\^ 
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Yi^l - 20.02 
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Ret* Coram. - Loc. 067 of Line 1 

If Xj^<X<Xj^^j_, then approximatfljr 3-f21c 
revolutions = OO-l-SSk ms« 



Error Stops? . . . . ^ . . If xoc^, routine will ring bell and halt. 
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^ 
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FORM 

BCD 

108-6.0 



SPECIFICA^riOIg 

lype:* . . ♦ Subroutine for G-l^D 

Mode of Operation:. •♦♦♦♦♦. Mxed point, sSinjle precisj.on 
B^xecution: Command Line 1 

Entry: 

(a) Multiplication At word time OUO 

(b) Divisj.on *At word tiine 035 

Exit: Normal or Overflow Ret. Comm. fi*om AE. at 

word time O38 

Scaling:. •. Input and output at decimal scale 10*^ 

Data Input: lO^^d >- 22.00 

10" Jc > 22.01 

10-2b > 22.02 

10-2a^ ->- 22.03 

Ret. Comm. for normal return to Line — ^21.02 
Ret. Comm. for return to Line after 
overflow ^22 . 03 

Data Output: lO^^y ^ 20.00 

10-2x = 20.01 

Normal Ret. Comm, » 21o02 
Overflow Ret. Comm. » 21.03 

Input data in Line 22 is preserved after 
multiplication only. 

Execution Time: 

(a) Multiplication II6 msec (U rev.) 

(b) Division • « .261 msec (9 revo ) 

Error Stops: None 

REMRKS 



Programmer must reset overflow before entering routine. Sj.nce Line 22 js 
used, magnetic tape write must not be :in process during execution of routine. 
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SPECIFICATIOHS 

^05* ..*..••.♦ Subroutine for G-1!>'D 

Mode of Operationr* . . .Fixed point, single precision 

Execution? , . .From Command Line 1 

^^"^^^y* .First command location = 073 

^^i't^ • . .At \vDrd time 054 

Scaling^- x and a must bo given with equal scale 

factors. 

The exponent "in" is to be scaled at 
2-^, However, it can be changed to 
another scaling by changing the increment 
in 074 to the same scale factor, x^ is 
obtained at the scale factor with iriaich 
A^ is entered. 

Data Input: a^ stored in L21 words 000 and 001 

m stored in L21 word 002 

a stored in L21 word 003 
X stored in AR 

Return stored in LI, word 054 

Data Outputr. ..... .x^ stored in L21, word 001 

Execution Timer 3 drum revolutions per iteration. Total 

time requirement is dependent upon the 
degree of precision desired before 
termination of calculation, and the choice 
of input constants. 

Error Stops r None 

REMARKS 

The choice of the constant "a" is made with consideration for the ease 
in solving a set of numbers to the same exponent and the need for as few 
iterations as possible in getting the desired answer. The exponent "m" 
will usually be an exponent of nonintegral value. 
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INTERCOM SUBROUTINE Sine 
T I t L E 



Prepared by : William Krause 



LOS ANGELES 45. CALIFORNIA AppfOved fay '. BENDIX RADIO DIVISIOH 



Paie J^ol 2 
Date; 3-13-57 
No. 



APPROXIMATIONS USED: 



Function: Sin V/S x 



Range: 
Approx : 



-l^x ^1 



Sin 



7f/2 I = Ct^ X + C3 x3 + C^ x^ + C^ x"^ 



Cj^ 1.57079UB52 

C3 -.61^5920978 

C^ .079U87663 

Cy -.00U362U76 



Max. Relative Error: 



.000001 



Function: Arcsine x^ 77^2 - Yl^' M'(x) 



Range: 
Approx 1 



O^x ^1 



^«(x) - a© * ai X + ag 



x^ + ao x3 + a,, x*^ + 



a^ X 



^0 


1.570795207 


*1 


-.21U512362 


ag 


.087876311 


^3 


-.OUU95888U 


% 


.0193U9939 


*5 


-.OOU337769 



Use is made of ti^i,". ivU.;.-iui ■'(■'•;; ; 



FORM 

BCD 

fOI-6.0 



cos 7//2 X = sin 'ti/l (x + l) for cos routine 



INTEaCO!] SUiiitOUTii'JE uiue 



TITLE 



Date : 

3-13-57 



Paie ^ol 2 

No. 



APPROXH^ATIOIIS USEDt Continued 



Sin e 



tan 9 



"Xl 1 + tan2 e 



for arctan routine 



Greater accuracy of the routines may be achj.eved by inserting the 
constants via an interdom program prior to punching the routine on 
tape. For example, to obtain C^, use programs 



Clear and add 


51.15707 


Add 


li6.9i;852 


Store in Loc 


036 



If it is known that the angle (of which the sine or cosine is to be 
obtained) lies wi.thin a certain iSo"" range, the routine may be speeded 
\£ considerably by dividing by appropriate constant {^ /2 or 90 ; after 
putting angle in principal value range (O^cos r: 180 or 
-90 ^ sin e - 90), then entering at 015 for cos and 017 for sin. 



ENTRIES 



sin 9 Idegi'cei? 000 



cos 9 



arctan X 
arc sin X 



\ radius 001 

! degrees 003 
radj.us OOU 

037 



Inputs e or X to A regj.ster 

Output t X or 9 in A register 



sin ent « 0^1 
cos ent ^ 3iU 
arctan ent « 37 
arcsdji ent « U; 



Reference: Hastings Approx, for Digital Computers. 
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5v 
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5v 


76 









Multiply '\ 1 - X 
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Subtract 7r/2 


63 
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Clear and add x 
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Transfer if neg. to 69 
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Clear and Subtract - 
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li9 
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Retvirn Transfer 
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Ux 


00 









Clear and Add 9 
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Uv 


79 









Return I^ansfer 
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00 
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»3 
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ai 
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Multiply by x 
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3U 
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35 


- 


50,6U592 
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51.15708 
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5v 
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77 









Divide by A 1 + y^ 
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78 









Store X 
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77 









Take abs. value of x 


U5 


13 


00 









^Store 1 X 1 


U6 
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78 









Negate 


U7 


29 


00 









Add 1 


kS 


^v 
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Square Rt* 
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